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ABSTRACT 

X-ray luminosities and surface brightness profiles of the hot gas haloes of simulated 
disc galaxies at redshifts z — — 2 are presented. The galaxies are extracted from fully 
cosmological simulations and correspond in mass to the Milky Way. We find that the 
bolometric X-ray luminosities of the haloes decrease by a factor 4 — 10 from z ~ 1 to 
z ~ 0, reflecting the decrease in the rate at which hot halo gas cools out on to the 
disc. At all redshifts, most of the emission is found to originate within 10-15 kpc of 
the disc. When combined with models in which the evolution of disc X-ray luminosity 
is dominated by X-ray binaries, the predicted halo luminosities at z ~ 1 show good 
agreement with constraints from spiral galaxies in Chandra Deep Field data. There is 
an indication that haloes with a metal abundance of 0.3Z Q overpredict observed X- 
ray luminosities at z ~ 1, suggesting that halo metallicities are lower than this value. 
Prospects for direct detection of the haloes of Milky Way-sized galaxies with current 
and future X-ray instrumentation are discussed. It is found that XEUS should be able 
to single out the halo emission of highly inclined Milky Way-sized disc galaxies out 
to z « 0.3. For such galaxies in this redshift interval, we estimate a lower limit to the 
surface density of detectable haloes on the sky of ~ 10 deg -2 . More generally, owing 
to their luminosity evolution, the optimum redshifts at which to observe such haloes 
could be 0.5 < z < 1, depending on their assembly history. 

Key words: methods: Af-body simulations - cooling flows - galaxies: formation - 
galaxies: haloes - galaxies: spiral - X-rays: galaxies. 



1 INTRODUCTION 

Standard models of disc galaxy formation require spiral 
galaxies to be surrounded by large reservoirs of hot gas 
which should be emitting at X-ray wavelengths and from 
whic h gas should still be a ccreting on to the disc at present 
(e.g. White & Frcnk 1991). The detailed properties of these 
haloes and their role in the formation and evolution of 
galactic discs remain largely unknown from an observational 
viewpoint, as the haloes have so far es caped direct detec- 
tion. Possible exceptions are NGC891 fereeman fc Houckl 
Il997l: IStrickland et all l2003al) and the Milky Way itself 
(e.g. lSidher et al.ll996l : lPietz et alll998ft . if neglecting cases 
where (i) the galactic discs show evidence of being disturbed 
by tidal interactions, or (ii) the X-ray emission at off-disc 
distances of a few kpc can be attributed to processes orig- 
inating in the disc such as feedback from star formation 
JPahlem et alJl2003t Strickland et al. 2003a,b). A possible 
explanation for this lack of halo detections may be pro- 
vided by recent simulation work feirnboim fc Deke]|l2003t 
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iKatz et aljl2003fi . which suggests that most galaxies accrete 
the majority of their gas at temperatures much lower than 
the virial temperature of their halo. This would be particu- 
larly pronounced at higher redshifts, and would imply that 
most of the halo radiation is emitted as Lya emission close 
to the disc rather than as soft X-rays from an extended, 
quasi-spherical region. 

In a previous paper dToft et alJl2002l . hereafter Paper I) 
we extracted disc galaxies from fully cosmological hydrody- 
namical simulations and determined their present-day halo 
X-ray luminosities, finding consistency with observational 
upper limits on halo emission from nearby spirals. It should 
be mentioned that we have subsequently detected an error 
in the X-ray calculations employed in that study. Correcting 
this increases the luminosity of hot haloes with gas metal- 
licity Z = 0.0 by a factor ~ 2 and that of Z ~ 0.3Z Q 
haloes by a factor ~ 1.25. This correction applies equally 
to all our simulated galaxies, so the overall conclusions of 
Paper I remain unchanged. Accounting for this error still 
produces results in consistency with observations, as the up- 
per limits on observed luminosities can easily accommodate 
the correction. One main result of the study of Paper I was 
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the conjecture that spiral haloes were possibly one order of 
magnitude brighter in soft X-rays at z = 1 than today, a 
result which was based solely on the predicted mass accre- 
tion rates of the discs. Analytical models of halo emission 
overpredi ct emission at z = by at least an order of magni- 
tude (see iBenson et aljl2000l and Paper I) and thus cannot 
be assumed to provide a reliable description of the evolution 
in X-ray emission with redshift. To amend this situation, we 
here extend our previous work by studying the high-redshift 
properties of a few disc galaxies extracted from cosmological 
cold dark matter simulations. 

We found in Paper I that the halo luminosity approx- 
imately scales as Lx oc V c 5 , where V c is the characteristic 
circular speed in the disc at R2.2 ~ 2.2Rd and Rd is the disc 
scalelength. Although subject to considerable scatter, such 
a trend is consistent with expectations from simple cooling 
flow models and suggests that massive galaxies (V c > 300 km 
s _1 ) would be optimum targets for observing halo emission. 
Such galaxies are rare, however, and hence do not compare 
well to the typical spiral seen in deep X-ray surveys. The aim 
of this work is therefore to study the predicted properties 
of, and the detection prospects for, haloes of more typical 
disc galaxies of size similar to the Milky Way (MW). 

The simulations and X-ray computations are described 
m S0 Results are presented in and compared to obser- 
vations in 2] We discuss the possibility for detecting the 
haloes in ^and summarize our findings in &\ 



2 SIMULATIONS AND X-RAY 
CALCULATIONS 

The primary goal is to study the redshift evolution in halo 
X-ray properties of simulated disc galaxies. To this end, 
we have extracted two disc galaxies from fully cosmologi- 
cal TreeSPH simulations of galaxy formation in a ACDM 
cosmology with Qm = 0.3, = 0.7, and Ho = 65 km 
s _1 Mpc -1 . These values are adopted throughout this pa- 
per. The simulations include star formation, stellar feedback 
processes, and a meta-galactic UV radiation field. For a de- 
sc ription of initial conditions and simulation details, we refer 
to lSommer-Larsen. Gotz fc Portinaril l)2003h . 

From a simulation with an adopted universal baryon 
fraction of ft = 0.10, we study in detail two galaxies (here- 
after denoted gall5 and gall8), each run both with a cool- 
ing function based on a primordial gas composition ([Fe/H] = 
—00, i.e. Z — 0.0) and with one based on the typical in- 
tracluster gas metallicity ([Fe/H]= —0.5, i.e. Z ~ O.3Z0). 
The latter abundance can probably a priori be considered 
a reasonable upper limit to the metallicity of hot halo gas. 
The two galaxies were selected from the requirement that 
they should show disc circular speeds V c at 2 = com- 
parable to that of the present-day MW. With values at 
z = of V c = 225 (250) km s" 1 for gal 15 with Z = 0.0 
[Z = O.3Z ) and V c = 202 (213) km s" 1 for gal 18 with 
Z — 0.0 (Z — O.3Z0), they have MW-like masses and hence 
correspond to L ~ L* galaxies in the local Universe. The 
z — version of both galaxies were included in Paper I. 
Since the galaxies were formed in the same cosmological sim- 
ulation, differences between them only reflect their different 
assembly histories. 

For the calculations of halo X-ray properties we follow 



the approach of Paper I, to which we refer for more de- 
tails: All SPH gas particles outside a box sized (1000 kpc) 3 
centred on the galaxy are cut away, as is 'cold' gas (with 
T < 3 x 10 4 K). Spatially smoothed temperature and den- 
sity fields are constructed from the values associated with 
individual particles, and X-ray volume emissivities are cal- 
culated from these using the meka p lasma emissivity code 
iMewe. Lemen fc van den Oordl ir986). 

At any given redshift, some variation in total X-ray lu- 
minosity is seen on time-scales of a few tens of Myr, due to 
galaxy-galaxy interactions and merging of smaller satellites. 
Since we are mainly interested in the long-term evolution of 
Lx, on the order of 10 frames separated by time intervals 
At ~ 100 Myr have been extracted around each redshift. 
This time interval is large enough that individual frames 
can be considered reasonably independent, given the typical 
time-scale involved, yet small enough for overlaps between 
adjacent redshift bins to be insignificant. X-ray properties 
of each frame have been calculated as prescribed above, and 
frames have then been omitted for which the total bolomet- 
ric (0.012-12.4 keV) luminosity deviates by more than 2a 
from the mean at that redshift (discarding on average less 
than one frame per redshift, corresponding to ~ 7 per cent of 
all frames). Mean values and dispersions of all quantities of 
interest were derived at a given redshift from the remaining 
frames and used in the subsequent analysis. 

For the Z = 0.0 galaxies we follow the halo properties 
back to z = 2.9 (gall5) and z = 2.3 (gall8), beyond which 
the disc itself is not well-defined. The Z = O.3Z0 galaxies 
show a more complex gas distribution at all redshifts and 
are followed only back to z ~ 2. 

3 RESULTS 

3.1 Lx versus redshift 

Of prime interest is the evolution with redshift of the X-ray 
luminosity, in part because this is expected to reflect the 
gas accretion history of the stellar disc as will be discussed 
below. 

Fig-0shows the redshift evolution of the rest- frame 0.2- 
2 keV and bolometric luminosities Lx,boi of the hot haloes of 
the two galaxies. In general Lx,toi is seen to increase with 
redshift, rising fairly steeply out to z ~ 1, beyond which 
the increase levels off. The increase from z = — 1 is a 
factor ~ 4 — 10, thus verifying our prediction from Paper I 
that disc galaxy haloes could be up to an order of magni- 
tude brighter at z = 1 than present-day haloes. As can be 
seen, gall5 shows a steady evolution in both X-ray bands, 
whereas gall8 exhibits a slightly more complex behaviour. 
The reason is that gall8 experiences a period of enhanced 
accretion, as will be demonstrated below. 

It is seen that at high z the 0.2-2 keV emission is some- 
what higher for the Z = 0.3Z© hot haloes than for the 
Z = 0.0 haloes, falling short of the bolometric luminosity 
at all 2 only by a factor of ~ 2 — 3. This is because the 
Z — 0.3Zq haloes display higher emission-weighted mean 
temperatures (~ 0.20-0.25 keV, with little redshift depen- 
dence) than the Z = 0.0 haloes (~ 0.10-0.15 keV), thus 
radiating a relatively larger fraction of their total emission 
in this band. Given the near-constancy of the hot gas tem- 
perature, the decrease in Lx with time is predominantly due 
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Figure 1. Redshift evolution of the rest-frame bolometric and 0.2-2 keV X-ray luminosities of the simulated galaxies. Top panel gall5, 
bottom panel gal!8, with metallicities as labelled. Error bars signify the la dispersion between frames at a given redshift. 



to a decline in the available amount of hot halo gas along 
with a decrease in the volume-weighted mean density of this 
gas. The mass of hot halo gas in the X-ray emitting volume 
declines only by a factor ~ 2 (Z = 0.0) and ~ 4 (Z = 0.3Z Q ) 
from z ~ 1 to z ~ 0, however. An implication of this com- 
bined with Fig. is that haloes radiate relatively more per 
unit mass of hot gas at high z than at z — (typically by a 
factor ~ 2 at z — 1 relative to z = 0). 

Note that for both galaxies, increasing the cooling effi- 
ciency by raising the gas metallicity to Z — 0.3Zq leads to a 
decrease in bolometric Lx at z = 0, consistent with the re- 
sults from our much larger sample in Paper I. The decrease 
is less significant for gall8, which, as will be argued below, 
is probably somewhat atypical of the galaxies studied in that 
paper. Also note that this decrease persists back to at least 
z ~ 2. This may at first seem counter-intuitive, given the 
higher X-ray emissivity of a Z = 0.3Zq plasma. The reason 
is that already at redshifts z>2a larger fraction of hot halo 
gas has cooled out to form stars in the Z = 0.3Z Q galaxies, 
leading at z « 2 to a factor ~ 2 deficit in the total amount 



and volume- weighted mean density of hot gas relative to the 
Z = 0.0 case. 



3.2 Accretion rates 

The halo luminosity can be linked to the rate at which hot 
halo gas cools out and accretes on to the disc, as illustrated 
by the following considerations. 

The cooling time is the gas energy density divided by 
the cooling rate, 



| NkT 
A^T' 



(1) 



where A is the cooling function, rin is the hydrogen number 
density, and N = (n e + ^2 th) is the total number density 
of electrons and all ion species i. The mass cool-out rate 
dM of gas with mass density p in a small volume element 
dV becomes 



dM < 



pdV 



2 



(2) 
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where p is the mean molecular weight, m v the proton mass, 
and equality corresponds to the case of absence of exter- 
nal heat sources and pdV work. Since the emission-weighted 
mean inverse temperature is 



/ An 2 H ±dV 



/Ar&dV 

and the bolo metric luminosity is Lx = J AnjjdV, we find 



(3) 



M < 
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In this simplified picture, we therefore expect M to (roughly) 
scale with Lx{l/T) ew . 

As was shown in Paper I, the gas in the simulations 
is essentially two-phased, so M is here estimated at any 
given redshift by considering the rate at which hot gas (T > 
2 x 10 6 K) cools out to the cold phase (T < 3 x 10 4 K), see 
ISommer-Larsen et al.l i2003l) . 

In Fig. El we have plotted M derived for the Z - 0.0 
galaxies along with the expectation from Eq. Q. There is 
a remarkably good agreement in shape, particularly for the 
steadily accreting galaxy gall5. Although the agreement is 
slightly less convincing for gall8, it is, however, clear that 
this galaxy undergoes a period of strongly enhanced accre- 
tion around z ~ 1.3 which affects its X-ray luminosity ac- 
cordingly. In both cases the accretion rate at z ~ 1 is found 
to be ~ 5 times larger than at present, decreasing e.g. for 
gall5 from 2.5 to 0.5 Mq yr -1 from z = 1 to z = 0. Based 
on Fig. El (left), an approximate relation describing this evo- 
lution in the redshift range z = — 1 is 



IogM(z) » 0.6z - 0.25, 



(•») 



with M in M© yr~ . The observed difference in normaliza- 
tion between this relation and Eq. Q is due to the simpli- 
fying assumptions underlying the latter. We speculate that 
the main reasons for the discrepancy are in part the neglect 
of pdV work and in part a geometrical effect related to the 
fact that cooling proceeds in a cooling flow, so considering 
a fixed volume for the estimation of M is a clear simpli- 
fication. The problem is highly complex, however, and we 
cannot exclude that other effects play a role too. 

Given the predicted evolution of disc accretion rates and 
halo bolometric X - ray lu minosities, it is interesting to note 
that lBarger et all feOOll) found that accretion on to super- 
massive black holes, as measured by the bolometric X-ray 
luminosity of active galactic nuclei (AGN), follows a very 
similar trend to that seen in Figs. and El The logarithm 
of the black hole mass accretion rate was found to grow 
roughly linearly with redshift out to z ~ 1-1.5, above which 
the increase levels off. It seems reasonable to assume that a 
positive correlation should exist between the disc accretion 
rate and that of any central black hole (given the observed 
correlatio n between black hole mas s and host galaxy bulge 
mass, e.g. iMcLure fe Dunlodl2002D . in which case the ob- 
servations oT lBarger^t^lTl200^ir could be taken to support 
the overall trends predicted in Figs. and El 



4 COMPARISON WITH OBSERVATIONS 

Although only a few relevant observational constraints on 
halo emission are available even in the nearby Universe, we 
showed in Paper I that Lx of our simulated galaxies at z — 
agree with observed upper limits for haloes of nearby mas- 
sive spirals and with estimates of the amount and luminosity 
of hot halo gas in the Milky Way (these conclusions all re- 
main valid after correction for the error mentioned in iQ. 
Some recent observational progr ess in this context include 
the results of iKuntz et alJ 12003^1 on emission in the nearby 
spiral M101. These authors find that our halo prediction 
from Paper I for a galaxy like M101 (V c ~ 170 km s" 1 ) can 
easily be accommodated within the amount of diffuse emis- 
sion obtained for this galaxy, provided the halo gas density 
displays a sufficiently small v ertical scale he i ght. O ther re- 
cent results comprise th ose of iDahlem et all J2003I). St rick- 
land et al. (2003a,b), and lWang. Chaves fc irwinH2003l) . but 
in these cases the X-ray emission is either affected by tidal 
interaction with another galaxy or can be attributed to pro- 
cesses in the disc, and they cannot be directly compared to 
our models. 

Although at high redshifts there is currently no possi- 
bility of comparing the predicted levels of halo emission to 
direct observations, useful constraints can nevertheless be 
extracted from the Chandra Deep Field observations (these 
are superior in sensitivity to any other pointed X-ray ob- 
servation). Spiral galaxies, classified morphologically from 
Hubble Deep Field data, have been detected in the 1 Ms 
Chandra Deep Field-North (CDF-N) data down to a nomi- 
nal limiting 0.5-2 keV point source flux of ~ 3 x 10 -17 erg 
cm -2 s _1 , and using a stacking technique detection has bee n 
performed to even lower fluxes llHornschemeier et al.l bo02). 
The latter authors considered the integrated X-ray emission 
from L ~ L* Sa-Sc spirals at z — 0.4 — 1.5 which did not 
harbour a bright AGN or display strong starburst activity. 
Resulting rest-frame 0.5-2 keV luminosities in several red- 
shift ranges were derived and can be taken here as observa- 
tional upper limits to our predicted halo emission. Account- 
ing for disc emission, not included in our simulated galaxies, 
strengthens the CDF constraints on halo emissio n further. 
To this end, we consider the MW disc, for which IWarwickl 
( 2002) lists a total 0.5-2 keV luminosity of ~ 5 x 10 39 erg 
s~ , including all stellar and diffuse disc sources. This pro- 
vides a measure of the contribution to be added to that of 
our simulated haloes at z — 0, in order to obtain an es- 
timate of the integrated X-ray luminosity of the simulated 
galaxies. However, the MW disc is dominated by X-ray bi- 
naries (XRBs) at a level of - 3 x 10 39 erg s" 1 JWarwickl 



2002). This contribution to the disc emission is unlikely to 
remain constant with redshift, as it is expected to roughly 
scale with the disc star formation rate (SFR). Globally, the 
latter is known to rise a factor of ~ 10 between z ~ 0-1 (e.g. 
iMadau. Pozzetti fc Dickinsonlll998l). a result wh ich also ap- 
plies to the CDF-N in particular lCohenll2003h . Based on 
an assumed such coupling between disc S FR and Lx of the 
XRB population, iGhosh fc White! (1200 ll) present a model 
for the redshift evolution of Lx from XRBs. Assuming a 
Madau-type SFR profile, these authors predict a factor of 
~ 5 (3) increase in Lx from XRBs at z = 1 (z = 2) relative 
to z = 0. We have adopted the expect ed XRB c ontribution 
to the integrated disc emission from IWarwickl i2002T) and 
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Figure 2. Total disc accretion rate M (solid line) of gas, : 

evolve d it with redshift according to the iGhosh fc White! 
( 2001) model. An addition from other disc and bulge X-ray 
sources of 2 x 10 39 erg s _1 , assumed constant for simplic- 
ity, is included in order to match the observed total MW 
disc luminosity. The resulting disc luminosity of an MW- 
like galaxy is compared in Fig. [3] (left) to the nominal limits 
in different redshift bins from the CDF-N sub-samples of 
galaxies with spectroscopically confirmed and purely photo- 
metric redshifts, respectively. Typical uncertainties on the 
observed X-ray luminosities are displayed for comparison; 
the lcr statistical errors are ~20 per cent (A. Hornschemeier, 
priv. comm.), very similar at the relevant redshifts to those 
of our si mulated galaxies. Adopting the mea n Lx in the 
sample of Sha plev. Fabbiano fc Eskrided i200ll) as represen- 
tative of the typical spiral Lx at z = (compa r able t o the 
MW value adopted here) . IHornschemeier et alJ (120021) note 
that the observed increase in Lx fr om 2 = to 2 = 0.6 is 
consistent with the prediction of the lGhosh fc White! <l200ll) 
model. When taking the CDF values as upper limits, Fig.|3] 
suggests that this consistency can be extended to redshifts 
z ~ 1.5, provided that Lx of the other disc X-ray emit- 
ting components, such as hot diffuse gas in the disc and 
bulge, does not increase dramatically between z w 0.5 and 
~ 1.5. Some increase in diffuse disc luminosity is probably 
expected, however, due to supernova driven feedback asso- 
ciated with the increase in star formation with redshift. 

ICohenl J20051) has recently shown that the CDF-N spi- 
rals out to z = 1.4 follow an Ljf-SFR relationship com- 
parable to that seen for local galaxies, albeit with larger 
scatter. This result s upports the approach employed by 
IGhosh fc White! J200ll) in their XRB model. Motivated by 
this, and by Fig. |3] (left), we compare in Fig. |3] (right) the 
CDF c onstraints to our halo p redictions, having added the 
scaled IGhosh fc White! <l200ll) model (plotted as upwards 
arrows) to the upper limit of the highest predicted halo lu- 
minosity at any given redshift. Uncertainties on the CDF 
values are here omitted for clarity. As can be seen, the 
Z — 0.0 galaxies are easily consistent with the CDF-N con- 
straints, whereas the metal-rich galaxies are much closer 
to the observed values. When adding the IGhosh fc White! 
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the quantity Lx bol(^-/T), normalized according to Eq. I'll . 

teOOll) model to the predicted halo luminosities, the agree- 
ment with the CDF-N results is fairly good, our predictions 
for the integrated emission then being consistent with results 
from the photometric CDF-N sample. 

There is an indication, however, that the Z — 0.3Zq 
haloes, in particular gall8, could be slightly too X-ray lu- 
minous at the relevant redshifts. While this is certainly not 
a large excess, we see at least two possible explanations for 
this. Possibly gall8 is peculiar. It certainly displays a high 
Lx at z = for its circular velocity and metallicity when 
compared to other galaxies of the study presented in Paper 
I, supporting this idea. Alternatively, Fig. |H] may suggest 
that assuming halo metallicities of Z = O.3Z0 at all red- 
shifts probed here is unrealistic, overpredicting the result- 
ing X-ray emission at the CDF-N redshifts. The fact that 
gal 15, which at z — is very typical of the galaxies studied 
in Paper I, is also close to the CDF constraints at z > 1 
would seem to substantiate this possibility. While only in- 
dic ative, such a res ult would not be surprising; as argued by 
e.g. lRenzinil (120031) . the global metallicity of the Universe at 
z ~ 3 is likely to be ~ O.IZ0 rather than ~ 0.3Z©. On the 
other hand, for haloes at z = 1, the predicted unabsorbed 
0.5-2 keV flux originating inside a region lxl arcsec 2 on 
the sky (within the central ~ 10 x 10 kpc of the halo) is 
~2x 10~ 19 and ~ 6 x 10~ 19 erg cm~ 2 s" 1 for the Z = 0.3Z Q 
haloes of gall5 and gall8, respectively. This is well below 
the CDF-N point source flux for norm al spirals at this red- 
shift, ~ (3 — 6) x 10 -18 erg cm -2 s _1 (IHornschemeier et alJ 
2002). Note, however, that we have not attempted to quan- 
tify any systematic uncertainties on the predicted integrated 
emission (related to the choice of SFR model, the coupling 
between Lx and SFR, the assumed XRB fraction of disc 
Lx etc.). For example, it is not obviously appropriate to 
use a Madau-type SFR profile to extrapolate the SFR of 
individual galaxies to high redshifts. Taking again the MW 
as an example and neglecting variations caused by short- 
term starbursts, evidence suggests that the MW disc has 
experienced a roughly co nstant SFR back to at least z ~ 1 
jRocha-Pinto et all [2 000). Assuming a constant SFR, the 
IGhosh fc White! l2CKHl) ""model predicts a steady decline in 
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Figure 3. Left: Chandra Deep Field-North results (0.5-2 keV) on L fs L* gala xies, along with the predicted evolut ion in disc X-ray 
luminosity of a MW-like galaxy, assuming that X-ray binaries drive the evolution (Ghosh & White 2001; Warwick 2002). Solid line is the 
CDF-N subsample of galaxies with spectroscopic redshifts, dashed line the subsample with photometric redshifts only, with typical la 
errors overplotted. Right: CDF-N results and our predicted intrinsic 0.5-2 keV halo luminosities. Dotted line is gall5, while dot-dashed 
line is gall8, with metallicities as labelled. Upward arrows show the shift in Lx resulting from adding the disc luminosities of the left 
figure. 



disc Lx from z = to z = 2 of ~ 30 per cent, in which case 
our predictions for the Z = O.3Z0 galaxies do not violate the 
CDF constraints, but, in fact, show even better agreement 
with these values. 

W e finally note that when adding the lGhosh fc White! 
( 2001) model to the Z = 0.0 haloes, the integrated emission 
is well below that of the CDF values and so is easily consis- 
tent with these when viewed as upper limits. Predicted and 
observ ed luminosities can in this case be reconciled if us- 
ing the lGhosh fc White! (l200lT) model along wit h a different 
SFR p rescription, e.g. the 'hierarchical' model of lBlain et alJ 
(1999). 



5 PROSPECTS FOR HALO DETECTION 

In order to assess the possibility of detecting the X-ray 
haloes of MW-sized discs with current and future X-ray in- 
strumentation, we have computed X-ray brightness profiles 
perpendicular to the discs of the simulated galaxies. The 
discs have been oriented edge-on, and observer-frame sur- 
face brightness (energy flux per unit solid angle) and flux 
profiles as a function of off-disc distance have been calcu- 
lated in 40 kpc wide and 5-8 kpc high slices parallel to the 
disc. Spectra were calculated for each SPH particle in the 
rest-frame band matching the observer-frame band of in- 
terest, and the resulting total surface brightness has been 
diminished by a factor (1 + z) 4 to account for cosmological 
dimming. Galactic absorption has been incorporated assum- 
ing an absorbing column density at z = of Nh = 3x 10 20 
cm" 2 , using the photoel e ctric absorption cross sections of 
iMorrison fc McCammor] jl983l) . Profiles have been calcu- 
lated in the 0.3-2 keV band, expected to provide the op- 
timum signal-to-noise (S/N) ratio for Chandra and XMM- 
Newton, given the predicted halo temperatures, as well as 



telescope effective areas and background levels. Calculations 
were repeated for the 0.1-3 keV band, which, in this con- 
text, is likely to be relevant to the next-generation XEUS 
mission, being well-calibrated down to energies at least as 
low as 0.1 keV and at high energies providing a much higher 
S/N than e.g. XMM-Newton. 

The 0.3-2 keV profiles are shown in Fig.|l]for the most 
optimistic case of Z — O.3Z0 edge-on galaxies. An immedi- 
ate conclusion, apart from the apparent fact that resulting 
fluxes are exceedingly low, is that at all redshifts the major- 
ity of emission originates within the innermost 10-15 kpc, 
as was also found to be the case at z = (Paper I). At least 
two additional inferences from these plots are worth noting: 

(1) Due to the strong luminosity evolution of the Z = 
O.3Z0 haloes (Fig. 0, the optimum redshift for observing 
such haloes in terms of their surface brightness is not nec- 
essarily z = but is not well-defined either. If Lx depends 
as strongly on redshift as for gall5, the optimum redshifts 
could be z > 0.5, depending on the assembly history of the 
galaxy. For the Z = 0.0 galaxies (not shown) there is a larger 
span between the flux at different redshifts, with z — al- 
ways being the optimum redshift. The fluxes for haloes at 
z > 1 are in all cases lower than those shown (this is also 
the case in the 0.1-3 keV band). 

(2) To unambiguously detect the haloes in imaging ob- 
servations, we estimate that a firm detection at least ~ 10 
kpc above the disc is required in order to avoid confusion 
with disc emission. Taking gall5 as representative, it can 
be inferred from the bottom panel of Fig. 2] that detection 
of haloes of MW-like galaxies at cosmological distances is 
clearly beyond the capabilities of current instrumentation, 
despite the expected redshift evolution of Lx ■ Even with the 
sub-arcsec spatial resolution of Chandra and point source 
sensitivities exceeding those of the Chandra Deep Fields by 
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Figure 4. Top panel: 0.3—2 keV surface brightness of the Z = 0.3Zq galaxies inside 40 kpc wide slabs oriented along the disc, as a 
function of vertical distance \z\ to the disc. The galaxies are viewed edge-on; for each vertical bin, the values above and below the disc 
have been added. Bottom panel: Corresponding flux profiles (surface brightness integrated over the sky area covered by the slabs), as a 
function of observed angular separation from the disc. 



a factor of a few, a direct detection of halo emission would 
not be possible. 

With an expected spatial resolution comparable to 
XM M-Newton, the next-generati on mission Constellation-X 
(e.g. IWhite fc Tananbaundll99g|) will not be able to image 
the haloes directly either. However, normal spiral discs in 
the nearby Universe typically exhibit two distinct spectral 
features, aT~ 0.2-0.3 keV thermal component from diffuse 
disc emissi on and a power-law from the point so urce popu- 
lation (e.g. lKuntz et al.ll2003l : ISwartz et al.ll2003ft . One may 
therefore rather seek to test the existence of the haloes and 
their predicted luminosity evolution by finding evidence for 
an additional low-temperature thermal component in the in- 
tegrated spectra of disc galaxies. To explore this possibility, 
we generated artificial source spectra for a Constellation-X 
observation of gall5 (z = 0.3), using xspec vll.O along 
with the response matrix for the calorimeter detector 1 . The 



Available from http://constellation.gsfc.nasa.gov/docs/ 



adopted response assumes a four-spacecraft configuration 
with a total effective area of 15,000 cm 2 at E = 1.25 keV. 
The source was modelled as a combination of a T = 0.3 
keV thermal mekal plasma with Z — Zq (representing disc 
diffuse emission), a T = 1.9 power-law (representing point 
sources), and a mekal plasma representing the hot halo, 
with each component normalized according to Figs. and 
[3] We find that a 1 Ms observation will gather around 300 
source counts over the full detector band (0.25-10 keV), of 
which only « 10 will originate in the halo. These numbers 
will be down by an order of magnitude at z = 1. It thus 
seems impossible that Constellation-X should be able to de- 
tect the haloes within reasonable exposure times. 

From the 0.1-3 keV profiles (not shown), one finds that 
extended energy coverage does result in larger observed flux, 
but the increase in emitted flux when including lower ener- 
gies is nearly balanced by Galactic absorption. Depending 
on redshift, the increase in observed flux from 0.3-2 keV 
to 0.1-3 keV is at most a factor of 1.5-2 for the Z — 0.0 
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galaxies and even less for Z = O.3Z0, even for a relatively 
low column density Nh as the one adopted. However, the 
Wide-Field Imager onboard the final configuration of XEUS, 
with a spatial resolution comparable to that of Chandra and 
an improved sensitivity particularly at low energies (an ex- 
pected limiting point-sour ce flux of 4 x 10~ 18 erg c m" 2 s _1 
for a 100 ks observation; iBleeker fc Mended boo^ . should 
observe ~ 6000 source counts in 1 Ms for gall5 at z = 0.3, 
of which ~ 300 would be from the halo. While this is not suf- 
ficient to warrant the introduction of an additional thermal 
component in a spectral fit, the halo should be detectable as 
extended emission surrounding the disc in soft-band images, 
after removal of disc and obvious point source emission. We 
thus predict that XEUS should be able to single out halo 
emission of highly inclined disc galaxies out to z ~ 0.3, but 
it will require exposure times comparable to those of the 
CDF. At z fa 1, about 120 halo counts out of 300 in total are 
expected. Although the ratio of halo to disc emission is sig- 
nificantly higher at this redshift, the lower flux and smaller 
spatial extent will make halo detection more difficult. 

The use of a different prescription for the disc SFR 
should not impact on the halo detection prospects, but the 
inclination of the galaxy will. The above results all apply to 
edge-on galaxies, i.e. a galaxies in which the polar axis is 
inclined by i — 90° with respect to the line of sight. For a 
galaxy with smaller inclination, halo detection will be more 
difficult, the effect in general depending on the 3-D distribu- 
tion of hot halo gas. One may estimate the minimum galaxy 
inclination at which a halo is still detectable. As a first ap- 
proximation, we can simply evaluate the sky area covered 
by the disc and assume that this area will effectively block 
all halo emission. This seems a reasonable assumption, since 
all halo emission originating in front of the disc will be hard 
to distinguish from disc emission, while halo emission be- 
hind the disc will be subject to the same effect as well as 
suffer from severe absorption by neutral gas in the disc. Al- 
though the method only takes into account the area blocked 
by the disc, not the fact that the halo itself will show a sim- 
ilar inclination, the error implied should be negligible for 
large inclinations. Assuming a circular disc of radius 15 kpc 
and employing the Z — O.3Z0 halo of gall5 at different 
redshifts, we estimate that at an inclination of 80° around 
half of the halo emission will still be detectable, whereas at 
i = 75° this number is down to ~ 30 per cent. Adopting a 
tolerance level of 50 per cent, we thus require the inclination 
to be i > 80°, i.e. the disc should deviate less than 10° from 
an edge-on appearance for the halo to remain detectable. 
This implies that ~ 10 per cent of all MW-like discs at a 
typical distance of z = 0.3 should display a detectable halo 
in a 1 Ms XEUS observation. 

The haloes of massive (V c > 300 km s _1 ) disc galaxies 
are predicted to have X-ray luminosity and surface bright- 
ness at z ~ an order of magnitude larger than the haloes 
studied here (cf. Paper I). If the halo luminosities of such 
galaxies follow a similar behaviour with redshift to that 
shown in Fig.0 and the galaxy shows an inclination i > 80°, 
these haloes should in fact be detectable to redshifts z ~ 1 
and less inclined haloes to well beyond z ~ 0.3. We can es- 
timate the total surface density of detectable haloes on the 
sky by assuming that 10 per cent of all haloes around discs 
with V c ^ 220 km s _1 can be detected to z = 0.3. Usingj 
the velocity function of spirals derived by iGonzalez et alJ 



(2000), one arrives at a density of w 10 haloes deg -2 in the 
adopted cosmology. Given our assumptions this would be a 
conservative lower limit, since some haloes beyond z = 0.3 
would also be detectable (those with V c well beyond 220 
km s -1 ), as would some at z < 0.3 (where our inclination 
requirement, i > 80°, is severely relaxed). 



6 SUMMARY AND CONCLUSIONS 

From the X-ray properties of the hot gas haloes of two Milky 
Way-like disc galaxies extracted from a cosmological simu- 
lation, run both with primordial and intracluster chemical 
composition, we find that halo X-ray luminosities increase 
roughly an order of magnitude from z = to z = 1, evolv- 
ing more gently at higher redshifts. There is good agreement 
between the redshift evolution of Lx(l/T) and that result- 
ing from rough estimates of the disc accretion rate, as sug- 
gested by simple analytical considerations. The logarithm of 
the disc accretion rate increases approximately linearly with 
redshift out to at least z « 1, from a typical value of ~ 0.5 
M yr" 1 at z = to ~ 3 M yr" 1 at z = 1. 

When added to a constant contribution from diffuse gas 
in the disc and bulge along with a model predicting the evo- 
lution of Lx from X-ray binaries, we find that the halo Lx 
of the galaxies is consistent with values derived for spirals 
in the CDF-N data. The luminosities of the Z = O.3Z 
haloes suggest, however, that a constant halo metallicity of 
O.3Z0 across the z ~ 0-2 redshift range studied here over- 
predicts the level of X-ray emission. Although this result 
is tentative rather than conclusive and depends on the as- 
sumed evolution in disc star formation rate, we speculate 
that halo metal abundances would be somewhat lower in 
the 2 = — 2 range than the typical metallicity of the in- 
tergalactic/intracluster medium (IGM) at z < 1. Possible 
explanations could include enriched supernova ejecta falling 
back on to the disc (galactic chimneys), enriched gas being 
expelled out of the galactic gravitational potential (possi- 
bly through supernova-driven superwinds during periods of 
strong star formation), or enriched gas being stripped from 
the haloes during galaxy-galaxy and galaxy-IGM interac- 
tions. 

We have assessed the possibility for detecting haloes 
of MW-like galaxies and shown that observations of such 
haloes at cosmological distances will have to await a future 
generation of X-ray instrumentation. We find that XEUS in 
its final configuration should be able to detect halo emission 
out to at least z ~ 0.3 in a 1 Ms observation, and that 
the surface density of haloes detectable in such observations 
should be > 10 deg -2 . In terms of surface brightness, the 
optimal redshifts for detecting haloes could in fact be 0.5 < 
z < 1.0, owing to their luminosity evolution. 

It is worth emphasizing that these conclusions are based 
on a particular numerical simulation with a certain set of 
input values; the aim was to assess the amount of asso- 
ciated X-ray emission in this specific case rather than to 
test the effects of varying the physical parameters entering 
in the simulation (some effort was devoted to this in Pa- 
per I). Given the observational difficulty in directly detect- 
ing the halo emission, the best way to test and improve on 
the predictions presented here in the immediate future may 
indeed be through additional simulations, perhaps coupled 
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with advances in our still incomplete understanding of the 
formation of spiral discs. Higher numerical resolution and 
yet more realistic input physics with improved estimates for 
the important physical parameters such as baryon fraction 
and gas metallicity should refine the predictions presented 
here. In particular, simulations of disc galaxy formation in 
which the gas chemical evolution is self-consistently included 
may shed new light on the amount and physical state of hot 
halo gas around normal spiral discs (Sommer-Larsen, Porti- 
nari & Romeo 2003, in prep.). With the observed integrated 
X-ray emission from spirals acting as a basis, improved mod- 
els for the redshift evolution of disc X-ray luminosity would 
alternatively allow tighter constraints to be imposed on halo 
emission. 
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